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Abstract 

The renormalization group applied to perturbation theory is ordinarily used 
to define the running coupling constant in the spacelike region. However, to 
describe processes with timelike momenta transfers, it is important to have a 
self-consistent determination of the running coupling constant in the timelike 
region. The technique called analytic perturbation theory (APT) allows a 
consistent determination of this running coupling constant. The results are 
found to disagree significantly with those obtained in the standard pertur- 
bative approach. Comparison between the standard approach and APT is 
carried out to two loops, and threshold matching in APT is applied in the 
timelike region. 
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I. INTRODUCTION 



A fundamental issue in quantum chromodynamics (QCD) is the behavior of the strong 
interaction running coupling constant as = gf/47r. The basic research tool is perturbation 
theory (PT) with its renormalization-group improvement |3J. In the QCD case in the limit of 
large momentum transfers Q, this approach provides a logarithmic decrease of the running 
coupling constant as ~ 1/ ln(Q 2 /A 2 ), where A is the QCD scale parameter that determines 
where the theory becomes asymptotically free. The study of the behavior of as outside of 
the asymptotic region is more difficult. It is known that the direct use of PT improved by 
the renormalization group leads to infrared instability of as and unphysical singularities, 
for instance, a ghost pole at Q 2 = A 2 . Unphysical singularities of a perturbative running 
coupling constant precludes a self-consistent determination of the effective coupling constant 
for timelike momentum transfers. Recently, a new method has been proposed @ for con- 
structing the QCD running coupling constant in such a way as to retain the correct analytic 
properties. This method is called analytic perturbation theory (APT). The main purpose 
of this paper is to analyze the region of timelike momentum transfers on the basis of APT 
P||, and compare the results of the PT and APT approaches. 

It is well known that a theoretical description of important timelike processes such as 
e + e~ annihilation into hadrons, or of decay widths of the r-lepton and Z-boson into hadrons, 
requires analytic continuation of the running coupling constant from the spacelike (Eu- 
clidean) region of momentum transfers (q 2 = —Q 2 < 0) into the timelike (physical) region 
(g 2 > 0). Although this problem has been studied since the 1970's (see, e.g., ||), it still 
remains a subject of great interest (see, e.g., It is obvious that information on the 

running coupling constant obtained from timelike processes, for instance, from e + e~ anni- 
hilation into hadrons, corresponds to knowledge of the coupling constant extracted from 
spacelike processes such as deep inelastic scattering, if the transition from the Euclidean 
into the physical region is performed in a correct manner (see HH) without violat ion of an- 
alytic properties of the hadronic correlator IT(g 2 ) and the Adler function D(q 2 ). When the 
analytic properties are not respected, the question arises: To what extent does this breaking 
of analyticity influence quantities extracted from physical processes? It is impossible to 
answer this question within the framework of standard perturbation theory. On the other 
hand, the APT method retains the correct analytic properties of the Adler D-function and, 
in addition, gives simple analytic expressions that can be compared with corresponding PT 
expressions and, therefore, allows quantitative analysis of the influence that the breaking of 
Q 2 -analyticity has on the running coupling constant. 

The organization of this paper is as follows: In the following section, we discuss the 
procedure of analytic continuation from the spacelike (Euclidean) to the timelike (s-channel) 
region. In Sec. Ill we examine this procedure in the conventional one-loop PT approach, 
and demonstrate that it is inconsistent. In Sec. IV we resolve this problem through the 
APT approach, and in Sec. V we compare the results of these two schemes. We move on 
to two-loops in Sec. VI, and demonstrate the stability of the APT approach. A matching 
procedure for timelike momentum transfers is given in Sec. VII, where we show how the 
coupling constant depends on the number of active flavors. A summary of our results is 
given in the Conclusions. 
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II. EFFECTIVE COUPLING CONSTANTS IN THE TIMELIKE AND 

SPACELIKE REGIONS 



First, we note that in the standard approach, the running coupling constant in QCD as a 
function of Q 2 is determined by the renormalization-group analysis in the region of spacelike 
momentum transfers. However, to parametrize many physical processes, one needs to know 
the coupling constant in the timelike region. To be specific, many experimentally measured 
ratios R a , where, e.g., a = e + e~, r , Z . . . , can be written in the form R a = Ry(l + A CT ). 
Here A CT is a QCD correction and represents the parton level of description of a given 
process with electroweak corrections. A CT can be expressed through the imaginary part of the 
hadronic correlator 7Z(s) = Im n(s)/7r. To parametrize 7Z{s) in terms of QCD parameters, 
a special procedure of analytic continuation is required. With that end in view, one usually 
employs a dispersion relation 

D(z ) = — z — : — = —z as k. 1 

y ! dz Jo ( s - z f K > 

where z = q 2 , and the inverse relation 

x 1 r s+it . du ( z ) 1 r s+ie . d(z) 

Z7T1 Js-ie dZ ZTtl Js-ie Z 

where the contour joins points s =F i e and lies in the region of analyticity of the function 
D(z), going around the cut Re z > 0. 

We define the effective coupling constants a eff and df 1 , respectively, in the spacelike (t- 
channel) and timelike (s-channel) regions, using the notation a = a /An and dimensionless 
(in units of the scaling parameter A) momentum variables, by 

D(z) oc 1 + dx a(z) + d 2 a 2 (z) + ■ ■ ■ = 1 + d x d cS (z) , (3) 
Tl(s) oc 1 + n d s (s) + r 2 d 2 (s) + • • • = 1 +n df(s) . (4) 

Relations (|T|) and (0) and the equality d\ = r\ result in the connection between these effective 
coupling constants 

I rs+ie rl 7 

sfW = "^/ (5) 

Z7T1 Js-ie Z 

d° s (z) = -z r-^afis). (6) 

•'O [s — ZJ 

Therefore, the QCD corrections for the class of physical processes considered with time- 
like momentum transfers are to be parametrized, according to Eq. @, by the effective 
coupling constant a^ ff (s), which is explicitly related to a eS (z) by Eqs. (Q) and (^). 

In any finite order of PT, the analytic properties of the running coupling d(z) should be 
the same as for the effective coupling constant d cS (z). Therefore, the connection between 
t- and s-channel running coupling constants, a(z) and d s (s) is defined by equations like @ 
and (^|). In the one-loop approximation the effective coupling constants coincide with the 
running coupling constants and in higher loops, the connection depends on the physical 
process. 
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III. PT ANALYSIS 



Consider the above procedure of analytic continuation within PT. In the one-loop ap- 
proximation, the running coupling constant is of the form 

a PT (^) = ~ * 2 . = ^r^—,, z = ~Q 2 /A\ (7) 
Po In (Q7A 2 ) fa In (-z) 

where @o — 11 — 2n//3 is the one- loop coefficient of the /5-function corresponding to riy active 
quarks. Inserting Eq. ([?]) into Eq. fl5|) we obtain the following expression for the running 
coupling constant in s-channel 

-pt/ \ 1 f s+le dz 1 1 / 7r lns\ , . 

a, (s) = — / — - — = — — h arctan . (8) 

s v ' 2mp Js-ie z ln(-z) nf3 \2 tt J K ' 

This expression is physically meaningless, because it is negative for any s and does not have 
the correct asymptotics, that is, going as 1/ Ins as s — > 00; the reason will be explained 
below. 

There is another way of calculating a s , based on the Shankar method ||. Using ana- 
lyticity of the Z)-function in the complex z-plane with the cut along the positive real axis, 
we may pass from the integral along the cut, given by expression @, that is, around the 
contour C\ (see Fig. |1|), to an integral around a circle of radius \z\ = s in the complex 
z-plane, contour C2, parametrized by z = — s exp(i^), —tt < ip < tt, to arrive the expression 

-"PT"/ \ 1 f dz _ pT 1 dip 1 7T 

a s \ s ) = / — a \ z ) = / — = —5- arctan - — ; (9) 



2ni Jc 2 z 2tiPq J-w In s + \p> 7r/3 In s 

this is positive when s > 1 and possesses the correct ultraviolet asymptotics. It is just this 
expression that is used as a one- loop PT result for all timelike momenta s > 0: 

a s PT (s) = — arctan- — . (10) 
Po In s 



It is obvious that Eq. ( |10"D provides the restriction a^' PT "(s) < 2tt//3 for any s (see [fLO] 



Thus, a formal conversion of the PT one-loop running coupling constant in the spacelike 
region (0) into expressions for the coupling constant in the timelike region leads to contra- 
dictory results (H) and @. The reason can easily be understood if one applies the Cauchy 
theorem (see Fig. 1) to establish the connection between the integrals in Eqs. (Q) and (^j), 

1 r 1 r 1 

— / dz^p{z) = -— dzip(z) + res [ip(z),-l] , ip(z) = —— -, (11) 

2m Jc 2 2 r K\Jc x zm[—z) 

which is consistent with Eq. (|) and @ because the residue of the function ip(z) at the point 
z — — 1 is 1. Therefore, the discrepancy between Eq. (|8|) and Eq. ([|) is due to an unphysical 
ghost pole in Eq. (0) at z = — 1 that violates the required analytic properties of the running 
coupling constant. The inclusion of multiloop corrections does not solve this problem but 
rather produces new unphysical singularities, as we will see in Sec. VI. Therefore, keeping to 
standard PT approximations that violate the necessary analytic properties of the running 
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Z-plane 




FIG. 1. Integration contours in the complex z-plane. 



coupling constant makes it impossible to pass into the timelike region in a self-consistent 
way. This can, for instance, be demonstrated by making an inverse transition from the 
timelike into the spacelike region with the help the dispersion relation (fl|). Substituting the 
running coupling constant d s (s) given by Eq. (|9|) into integral (H) following from Eq. ([[]) and 
taking account of the expression arctan(7r/lns) = sgn(lns) 7r/2 — arctan(lns/7r) , we arrive 
at the formula 



"PT 



1 



ln(-z) 



1 



1 



1 + 2 1 - Z 



(12) 



which is different from the starting point, Eq. ([?]). 



IV. APT ANALYSIS 

The problem of how to make the correct transition between the space- and timelike 
regions can be solved in the framework of the APT method 0,[3| that ensures the correct 
analytic properties of the coupling constant without introducing extra parameters. The 
resulting one-loop expression for the analytic coupling constant in the Euclidean region is 
as follows: 



-APT 



1 

Jo 



In (-z) 



l + z 



(13) 



The first term in brackets determines the asymptotic behavior at large momenta and is of 
the form given in PT. The second term, of a nonperturbative nature, compensates the ghost 
pole at z = — 1. When one employs the analytic coupling constant (jT3P, both methods of 
calculating a s considered above produce the same result, i.e. 



1 



dz 



APT 



27ri JCi z 



1 



dz 



APT 



27ri Jc 2 z 



1 / 7T In s \ ... 

^ = ~~R ~ arctan ~ ' 
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We note that APT and PT give distinct values for the running coupling constant in the 
timelike region, for example, 



-APT 



(15) 



where af T is given by Eq. ([8]). Consistency of the APT approach also follows from the 
fact that we can reconstruct the initial expression ([13]) when the timelike coupling (|14]) is 
substituted into Eq. ([]). It is of interest to note that this consistency is due to the second 
term in Eq. (13) that compensates the pole, whose contribution to the integral around the 
contour C*2 equals zero when s > 1, i.e., we have the equality 



dz 
— i 

c 2 z 



APT 



dz 
— < 

c 2 z 



PT, 



z), S > 1 , 



(16) 



where the function a (z) is defined by Eq. (^). Therefore the PT expression (|?p gives 
the same result as the APT approach in the timelike region for s > 1 if the contour C2 
is used. However, there is no inverse correspondence for PT [see formula (12)]. Moreover, 
note that an equality analogous to Eq. ([UJ) does not arise if the integrand contains the 
running coupling constant multiplied by a function of z. For the i? T -ratio, for instance, a is 
multiplied by a polynomial in z and, as is shown in |]llj, the contour integral over C2 in PT 
turns out to be different from that in the APT approach.^] 
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FIG. 2. The behavior of the running coupling constant calculated by different methods as a 
function of the dimensionless variable s = q 2 /A 2 . 



X A detailed comparison of inclusive r decay can be found in Ref. [12] 
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V. NUMERICAL COMPARISON OF RUNNING COUPLING CONSTANTS 



The results obtained are illustrated in a series of figures. (In these figures we take 
the number of quark flavors to be 3.) Firstly, we consider the region of small momentum 
transfers. Fig. 2 shows the behavior of the running coupling constant computed by different 
methods in the region —10 < s = q 2 /A 2 < 10. The solid line represents the APT coupling 
constant calculated by formula ( |13|) in the spacelike region (number 1 in the figure) and by 
formula (|T4] ) in the time-like region (curve 1*). Dots denote the "PT" coupling constant 
determined by Eq. (|10"D (curve 2*) and by Eq. ([12]) (curve 2). The dash-dotted line 3 
represents the PT coupling constant computed by formula (|7|) in the spacelike region, and 
curve 3* by formula (§) in the timelike region. (Incidentally, note that curves 2 and 2* 
vanish at the origin, which is beyond the resolution of this figure, while curves 1 and 1* 
approach the universal value 47r//3o = 1.40 at the origin.) 

As is seen from Fig. 2, the behavior of the APT running coupling constants (curves 1 
and 1*) is almost, but not quite, mirror-symmetric, and at s = the space- and timelike 
APT running coupling constants are both equal to the universal value 47r//3 (see @). The 
pairs of curves of the standard PT approach, 3 and 2*, or 3 and 3*, do not show analogous 
behavior of the running coupling constants. In the spacelike region the function a PT grows 
without limit (curve 3), whereas in the timelike region (curve 2*) it is limited to the value 
2n/ Curve 2 calculated with the coupling given by curve 2* in the dispersion integral 
does not reproduce the initial curve 3. 

Now consider the region where the running coupling constant a ~ 0.3, which ap- 
proximately corresponds to the mass scale of the r-lepton, M T = 1.78 GeV, a(M T ) = 
0.34 ± 0.04 [13]), defined in the spacelike region. It is known that the decay r — > hadrons 
is important for testing QCD, as it allows the most accurate determination of the running 
coupling constant at comparatively low energies (see the review ||14|| ). 

Fig. 3 shows the behavior of the running coupling constant versus the dimensionless 
variable s = \q 2 \/A 2 . Notation is the same as in Fig. 2; a dashed horizontal line corresponds 
to a = 0.34 . As it is seen from the Fig. 3, curves 1,2,3 that describe the spacelike region 
noticeably differ from each other. With increasing s, they begin, as they should, to approach 
each other, which is demonstrated on the top right of Fig. 3. Values of the parameter A 
calculated with the running coupling constants described by curves 1 , 2 and 3 are different. 
For example, the value of APT-function (curve 1), equal to 0.34 is achieved at sq = 41.5, 
which corresponds to A APT = 276 MeV. For PT-curve 3 s = 60.5 and A PT = 228 MeV. 
Note that for curve 2 the value a = 0.34 cannot be achieved at any value of s. For timelike 
momentum transfers, recall that curves 1 * and 2* as functions of the dimensionless variable s 
coincide when s > 1. However, they are characterized by different values of A, which results 
in different values of the running coupling constant in the timelike region, a^ PT (M r ) = 0.31 
and a'g PT " (M T ) = 0.29. (The procedure here is to use the spacelike value of a to determine 
Q 2 /A 2 , and then, with the same numerical value of s, determine a s from Eq. (|14]) and 
Eq. ([TI]), respectively.) With the accuracy attained at present for experimental data on the 
hadronic decay of the r [13]], this quantitative discrepancy is becoming significant. 

Let us finally observe the evolution of the running coupling constant in the region of 
momentum transfers of order of the Z-boson mass Mz = 91.2 GeV. The running coupling 
constant a, corresponding to the dimensionless variable s is drawn in Fig. 4. (At this point, 
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FIG. 3. Behavior of the running coupling constant as a function of the dimensionless variable 
s = \q 2 \/A 2 . Notation is the same as in Fig. 2; a dashed horizontal line corresponds to a = 0.34 . 
The graph on the right top shows the behavior ot the same curves for large values of s. 
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FIG. 4. Behavior of the running coupling constant in the vicinity of the Z-boson mass. Curves 
1,2,3 and 1* ,2* represent the same coupling constants as in Fig. 2 and Fig. 3. 
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we neglect the change in the number of active quarks with growing energy since this effect 
does not change the overall picture. We will consider the change in the number of active 
flavors at the two-loop level in Sec. VII.) The curve denoted by a represents all three curves 
1 , 2 and 3 that are drawn in Figs. 2 and 3 and describe the behavior of the running coupling 
constant in the spacelike region, which merge into one curve with high accuracy for these 
large values of s. The curve denoted by a s corresponds to the coupling constant in the 
timelike region and to curves 1* and 2* plotted in Figs. 2 and 3. In the region of sufficiently 
large values of s, the well-known approximate formula with the so-called 7r 2 -term (see, e.g., 

00) 



TV 2 



a, = a 



3 In 



2 



(17) 



works well (with an accuracy ~ 0.1%) both for PT and APT. This approximation gives a 
difference between a s and a of about 2%. Substituting the value of the parameter A fixed 
at q = M T , we obtain the corresponding values of the running coupling constant at q = Mz' 
a APT = 0.120, a PT = 0.117 (spacelike region); c^ PT = 0.118, a PT = 0.114 (timelike region). 
Thus, even at such large values of s, the effect of analyticity on the running coupling constant 
amounts to ~ 2%, i.e., it is comparable with the contribution from the 7r 2 -term and from 
higher PT loop corrections. 



VI. TWO-LOOP RESULTS 

We now extend the above considerations to the two-loop level. The distinction between 
the APT and the PT running coupling constants in the Euclidean region has to do with 
the unphysical singularities of the PT running coupling constant. Following the results of 
Ref . @ , we can write down the analytic running coupling constant in the form of a sum of 
the standard perturbative part and additional terms which compensate for the contributions 
of the unphysical singularities, a pole and a cut: 

a APT (z) = a PT (z) + Aa po i c (^) + Aa cut (z) . (18) 

For the two-loop perturbative running coupling constant, we use [|2|.|T2"f 

aPT W^ L + Bllll ( 1 + L/ Bl ) ' L ^<-^'4- < 19 > 

where Bi = (3i/(3q 2 , and (3i = 102 — 38rif/3 is the two-loop coefficient of the /3-function. 
Obviously, at large L Eq. (^) gives the standard PT expression as an expansion in inverse 
powers of L, 

a ^.-( 1 _ B -) +0 (J 3 ). m 

According to Eq. (0), the contribution coming from the unphysical pole is cancelled by 
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^ z) = W a TTl> (21) 

while the unphysical cut is removed by the following compensation term 



1 /•cxp(-Bi) d a n B 

Aa cut (z) = — — / 5 -. (22) 

ir/3 Jo a + z [lno- + B x ln(-l - lna/5i)] 2 + 7r 2 5 2 

To calculate the analytic running coupling constant in timelike region, we use the ex- 
pression in terms of the spectral density g(a) || 

af^(s) = - f° — g(a) , (23) 

71 Js (J 

where g(a) = Im a PT (— a— ie). The spectral density g plays a central role in the APT method; 
the spacelike running coupling constant, a APT , is also expressed through g as follows: 

d APT (z) = - r ^— g(a). (24) 

7i Jo a — z — ie 

As outlined above, independently of the order of approximation, the APT running cou- 
pling constants defined in the space- and timelike regions at Q 2 = and s = are both equal 
to the universal infrared limiting value 4n/ j3o |§, which is important to establish the stability 
in the region of the small momentum transfers. (This result is proved in [|HJ.) Consider the 
region in which the value of running coupling constant a ~ 0.35 (the r lepton scale). Fig. 5 
shows the behavior of the two-loop running coupling constant for the same interval of the 
dimensionless variable s as in Fig. 3. The solid line represents the APT coupling constant in 
the spacelike region, which can be computed from (|TBD; this is like 1 in Figs. 2 and 3. The 
dotted line denotes the APT coupling constant in the timelike region, computed from (p3|); 
this is like 1* in one-loop. The dashed line represents the spacelike PT coupling constant 
defined by formula (|T9D, which is like 3. The dash-dotted curve corresponds to timelike PT 
coupling constant constructed taking into account 7r 2 -terms, like Eq. (0), and analogous to 
2*. This figure demonstrates that, as in the one-loop case, there is a difference in behavior of 
all these constants. Moreover, the region in which the value of running coupling constant is 
about 0.34, is shifted to smaller s; the value of 0.34 for the APT running coupling constant 
is achieved at Sq = 8.6 which corresponds to A APT = 607 MeV and the value of the timelike 
running coupling constant a APT = 0.32. For the PT running coupling constant, Sq, = 18 and 
A PT = 419 MeV. Thus, while in the one-loop case APT parameter A, is 20% larger than PT 
value, in the two-loop case this discrepancy increases to 45%. 

In Fig. 6 we show the stability of APT result for the ratio of the space- and timelike 
APT running coupling constants. The solid line corresponds to three-loop case, the dashed 
line is two-loop, and the dot-dash is one-loop. (The one-loop result was already given in 
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2 In Fig. 5 we do not plot the three- loop result (see Q), because it practically coincides with the 
two-loop one. 
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FIG. 6. The ratio of spacelike and timelike values of the running coupling constant in the APT 
scheme, at one, two, and three loops. 
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Let us now to demonstrate that in the two-loop approximation, as in the one-loop case, 
the equality (|16D is valid as well. By using Eq. (|18D , we get 

2vriaf T = f -a APT (z) = f ^{z) + / -Aa polc (z) + f -Aa cut (z). (25) 



As follows from ([21]), the pole term has the same structure as in the one-loop approximation, 
and its contribution to contour integral equals zero when s > 1. One can also find that the 
contribution of the cut term (p2]) in Eq. (|25|) equals zero when s > exp(— As a result, 
the equality (0) holds for s > 1 as well, but it should be stressed (see above discussion in 
Sec. V) that this not mean that the APT and PT timelike coupling constants coincide with 
each other. 



VII. s-CHANNEL MATCHING 

Let us now discuss the issue of how the parameter A changes with energy as the number 



of active quark changes: A — > A n , (see [16] for further details.) The relationship between 
A„, and A n/+ i may be fixed by the matching conditions for coupling constant at "quark 
thresholds" (see, e.g., |T7[]) which is usually applied to the running coupling constant in the 
Euclidean region. The APT method opens the new possibility of performing the threshold 
matching in the physical region, where the number of active quarks can be associated with 
the energy threshold of quark pair production. It is important to note that any matching 
procedure of the coupling constant in the Euclidean region, for which one uses the condition 
of the type ReQ 2 > £ 2 M 2 (usually 1 < £ < 2), leads to a violation of analyticity of a(z). 
Although s-channel matching demands only continuity of the function d s (s) at the threshold, 
and not of its derivatives, the spacelike running coupling constant a(z) will be an analytic 
function of z, and due to the representation (|]) a(z) "knows", in principle, about all quark 
thresholds. Therefore, the APT method gives a more consistent definition of the running 
coupling constant and a natural way to perform the matching procedure. Within the APT 
approach, we will require that the timelike function a^ PT (s) should be a continuous function 
at the threshold points: 

af T [(Q nf+1 /A nf ) 2 ,n f ] = af T [(Q nf+1 /A nf+1 ) 2 , n f + 1] , (26) 



where Q nf is defined by the pole masses M q of quark pair. Taking into account Eqs. ( |T6"D 
and the results of the previous section that the unphysical singularities do not contribute at 
s > 1 to the contour integral we can rewrite Eq. ( P6| ) in the following form 



/ 



<lZ a PT (z,n f )= [ — a PT (z,n f + l). (27) 



z J z 

\z\ = (Q„ f +i/A nf ) 2 \z\=(Q nf+ i/A nf+ i) 2 

Therefore, the conventional matching condition 



3 For < nf < 6 we have 0.4 < exp(-^i) < 0.6. 
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<f T [{Q nf+1 /A nf ) 2 ,n f ] 



[{Qrif+l 

/A nf+1 ) 2 ,n f + l], 



(28) 



which one usually uses in perturbation theory, is modified and written down as the relation 
of the contour integrals, Eq. (|27|) . 

As an example, consider a change of the two-loop scale parameter A when passing through 
a quark pair threshold in PT and APT by using a(M T ) = 0.34 and the following values of 
pole c-, b- and t-quark masses M c = 1.6 GeV, Mb = 4.5 GeV, and M t = 174 GeV and 
Q 4 = 2M C , Q 5 = 2M b , Qe = 2M t . In the perturbative case, we find A PT = 419 MeV, A PT = 
338 MeV, A PT = 230 MeV, and A PT = 92.6 MeV, the ratios of which obey the well-known 



relations from Ref. |7J. In the APT case, we obtain A^ PT = 607 MeV, A^ PT = 471 MeV, 



A^ PT = 316 MeV, and A^ PT = 129 MeV. The ratios of these quantities are close to the 
perturbative relations; however, the PT and the APT values of A with the same number of 
active quarks differ by about 40%. 
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FIG. 7. QCD evolution of the running coupling constants (defined in the spacelike region) 
compared to experimental data. 



In Fig. 7, we plot results for the QCD evolution of a,s{Q), comparing the APT running 
coupling constant as discussed above which uses s-channel matching according to Eq. fl26|) 
with the standard PT running coupling constant [see Eq. (|20|) 1 which uses the matching 
procedure given by Eq. (p8|), starting from Mz down to Q = M T . Also shown on the graph 
is the experimental data measured by L3 collaboration ||18|| . As experimental input we use 
the average value as(Mz) = 0.1207 ± 0.0016 from Ref. [IS]. (In order to do not encumber 
the figure we do not plot the corridor of errors). 

We should give a little explanation of how the spacelike APT running coupling was 
calculated. First, we calculate A 5 from the measured value of a APT (M z ) = 0.1207. Then, we 
use this value of A 5 to determine qA pt for all s, and hence through the matching procedure, 
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determine A 6 , A 4 , A 3 . The spectral density that we used above at rif = 3 is determined for 
arbitrary rif by af PT through j^g(s, A n/ , rif) = —s dd s (s,A nf ,rif)/ds. An explicit formula 
for g as a function of s, A, and rif is given in Eq. (24) of Ref. ||12|| . Then, from the spectral 
representation (p4|) , we find the spacelike running coupling constant from 



f /-4A/2 rfo MM, 2 A a 

« APT (Q) = 4 / A 3 , 3) + / ; g(s, A 4 , 4) (29) 

JO S + Q z J AMI s + V 

/■ 4M t 2 c/s , , r 00 ds . . . 

+ / , g(s, A 5 ,5) + / g (a, A 6 ,6 

J4M 2 S + ./4M 2 S + 

It should be stressed that in Fig. 7 we have plotted the value of the QCD running 
coupling constant extracted by using the perturbative parametrization. However, this is not 
really self-consistent. For instance, in the case of the semileptonic decay of the r-lepton, 
to parametrize the process in the term of the QCD scale parameter A one usually uses the 
analytic properties of the running coupling, which are obviously broken by the perturbative 
approximation due to the unphysical singularities. Within the APT approach it is possible 
to maintain the required analytic properties and give a self-consistent description of the 
process ||12|1 . That, in principle, changes the value of the QCD running coupling extracted 
from the experimental data. Thus, the experimental points, plotted in Fig. 7 should be 
considered as illustrative only. Nevertheless, it is clear that if we use a normalization point 
with a large value of momentum, the curve of the running coupling constant corresponding 
to the APT method lies below than the corresponding PT line, which, from the point of view 
of the perturbative description, corresponds to a smaller value of A at low energy. The fact 
that low energy data prefer small values of the scale parameter A and that at the same time 
the high energy data prefer larger values of A has been emphasized in Ref. [|T^]. Thus, this 
apparent discrepancy may be understood in the framework of the APT method; however, we 
should mention again that one needs to perform a reanalysis of the low energy experimental 
data by using the APT parametrization, as in the case of r decay |]I2" |, in order to extract 
the QCD running coupling constant. 



VIII. CONCLUSIONS 

Let us briefly summarize our considerations. To determine the running coupling constant 
in the timelike region, we took advantage of APT because it provides a consistent procedure 
necessary for analytic continuation. It is to be noted that the APT method ensures not only 
correct analytic properties of the running coupling constant but also stability with respect 
to higher loop corrections, which is essential for the stability of our procedure of analytic 
continuation. This stability is provided, in part, by the universal infrared limit value of 
the running coupling constant at q 2 = —Q 2 — ► —0 that is invariant with respect to higher 
loop corrections. The proposed method of constructing the running coupling constant in 
the timelike region results in a function with the same universal infrared limit value when 
q 2 -> +0. 

Quantitatively, our analysis shows that the effect of analytic continuation can be asso- 
ciated with 7r 2 -terms only at very large momentum transfers of the order of the Z-boson 
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mass where the contribution of the 7r 2 -terms is small. At intermediate and, especially, at low 
momentum transfers it is important to take account of the correct analytic properties of as, 
which permits a consistent transition into the timelike region. The Q 2 -dependence of as is 
essentially different from the dependence of as in PT. Our analysis shows that the popular 
PT expressions for as as expansions in l/ln(Q 2 /A 2 ), containing nonphysical singularities, do 
not allow a self-consistent interpretation of information obtained from different experiments 
on the evolution of as outside of the asymptotic region. From our numerical estimates it 
follows that analyticity of the running coupling constant has great influence on the value 
of the parameter Aqcd extracted from experimental data and on the Q 2 -evolution of as- 
Note that these considerations are also important for the investigation of power corrections, 
which are now under intensive study (see, e.g., |20| ). The importance of power corrections 
in the APT scheme relative to perturbative terms naturally will be different than in the 
conventional approach. As in PT, the influence of quark thresholds results in a reduction of 
the scale parameter A as the number of active flavors nf increases. However, the importance 
of maintaining the correct analytic structure suggests that the required matching be made 
in the physical region. 

The APT method appears to be fruitful for studying the problem of analytic continuation 
of as into the timelike region. There is no doubt that extracting more detailed information 
from experimental data on timelike processes requires a more thorough theoretical analysis 
within APT including the estimation of the contributions from higher order processes, mass 
corrections, and so on. These will be considered in our subsequent papers. 
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